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Abstract. Recent measurement of a high millisecond pulsar mass (PSR J1614-2230, 1.97 ±0.04 Mq) 
compared with the low mass of PSR J0751+1807 (1.26 ± 0.14 Mq) indicates a large span of masses of 
recycled pulsars and suggests a broad range of neutron stars masses at birth. We aim at reconstructing 
the pre-accretion masses for these pulsars while taking into account interaction of the magnetic field 
with a thin accretion disk, magnetic field decay and relativistic 2D solutions for stellar configurations 
for a set of equations of state. We briefly discuss the evolutionary scenarios leading to the formation 
CO ■ of these neutron stars and study the influence of the equation of state. 
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1. Introduction 

The discovery of a 1.97 ±0.04 Mq rapidly-rotating neutron star (NS) ( Demorest et al. 2010t ) 



motivates the theorists to understand its origin and composition. P opulation of mil lisecond 
^ '. pulsars contains also a low-mass PSR J0751+1807 (1.26 ± 0.14 M g . lNice et al.ll2008h : a sta- 
1 tistical analysis of the measured and estim ated pulsar masses, in particular of the millisecond 
, ones, was presented in (jZhang et al.ll201ll ). Despite a striking difference in their masses, their 

spin periods and inferred magnetic fields are similar (B oc \J PP, see Table Q}. Both of them 
were spun-up by accretion from a binary companion. The evolution of such systems is typically 
composed of several stages: starting from the Zero Age Main Sequence (ZAMS), more massive 
■ primary (NS progenitor) enters the Red Giant (RG) phase and engulfs the secondary for a 
brief Common Envelope (CE) phase. Subsequently the primary, now as a helium star of a sub- 
qq ■ stantially lower mass, explodes as a core-collapse supernova (SN), forming the NS. Secondary 
star evolves from a ZAMS star into a RG and spins-up the NS via the Roche lobe overflow 
disk accretion. Depending on stellar masses, the resulting composition of the secondary, bi- 
^vq | nary orbital periods and semi-major axes may be very different, which in turn may result in 
7— I ■ additional stages of evolution: Intermediate Mass X-ray Binary (IMXB) stage, connected with 
the Terminal Age Main Sequence instabilities in the case of PSR J1914-2230, or secondary 
evaporation in a Low Mass X-ray Binary (LMXB) - so-called Black Widow scenario - in the 
^ I case of PSR J0751±1807 (see Tab le M for quantitati ve figures and |Pe Loore fc Doomlll992l 



Tauris fc van den Heuvelll2006l and lTauris et al.ll201ll for more details). 



2. Methods 

Hereunder we focus on the recycling/spin-up stage of the evolution leading to the formation 
of a millisecond pulsar. It is modeled by a sequence of stationary rotating configurations of 
increasing mass, calculated for a given equation of state (EOS). Rigidly rotating axisymmet- 
ric stellar configuration were obtained by means of the library for relativistic computations 
LDRENE. Increase of the total stellar angular momentum J is calculated by accounting for 
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Table 1. Selected parameters of millisecond pulsars PSR J1614-2230 and PSR J0751+1807 and their 
binary systems (orbital period Pb, eccentricity e and the comp anion mass Mwd 5 mass measurement 
errors reflect 1-a confidence levels). 



PSR 


M [Mq] 


P [ms] 


P [10~ 21 s] 


Pb [d] 


e 


Mwd [Mq] 


J1614-2230 


1.97 ±0.04 


3.15 


9.6 


8.7 d 


1.3 x 10" b 


0.500 ± 0.006 (CO WD) 


J0751+1807 


1.26 ±0.14 


3.48 


7.2 


6 h 


5.6 x 10~ 5 


0.12 ±0.02 (He WD) 



the transfer of the angular momentum from the thin accretion disk interacting with the NS 
magnctospherc: 

dJ K B 

where an accreted (baryon) mass dMj, carries a specific angular momentum Iq from the inner 
boundary of a disk; M is t he accretion rate and Kb deno tes a magnetic (braking) torque. We 
employ the prescription of iKluzniak fe Rappaport] ( 2007 ), improved by considering the rela- 
tivistic effects (the innermost stable circular orbit, Beiger et al. 2011al ) and using approximate , 
but accurate representation for orbital parameters of particles in the disk (jBeieer et al.ll2010h . 
Gradual magnetic field decay, which is c aused by the infalling accre ted matter, is taken into 
account as well; we adopt the relation of Shibazaki et al. (1980) and Taam fc van den Heuvel 



( 19861) that links the decrease of B with the amount of accreted matter only (actual evolution 
of th e polar magnetic field of an accreting NS can be more complicated, see lZhang fc Koiima 
20061 and references therein) 



Followin g EOSs were used: APR by lAkmal et al. | (I1998TI . DH bv lDouchin fc Haensell (|200ll) 
and BM bv lBednarek et al.l (l201l[ ). APR EOS (A18 + Sv + UIX* model) is a variational, non- 



relativistic many-body solution with rclativistic corrections (maximal mass of a non-rotating 
star: Mg£ = 2.19 Mq, R(M^) = 9.93 km). DH EOS is constructed using a non-rclativistic 
energy density functional based on the SLy4 effective nuclear interaction, designed to describe 
both crust and core in an unified way (M%£ = 2.05 M , R(Mg£.) = 10.0 km). BM EOS 
is based on relativistic mean field model, with a non-linear Lagrangian that includes quartic 
terms in meson fields, and two additional hidden-strangeness mesons a* and </> that couple to 
hyperons only (M%£ = 2.03 M Q , R(M^) = 10.7 km). 



Detailed description of methods and tests can be found in lBeiger et al.l (|2011al ). 



3. Results 

The recycling process for PSR J1614-2230 and J0751±1807 is presented in Fig. Q] and [2j 
respectively, and visualized in terms of accretion time r, the amount of accreted mass M acc , 
and stellar spin frequency /. We assume constant accretion rate M during the recycling. The 
lifetime of an LMXB or IMXB system puts strong limits on the a mount of accr e ted m ass and 
M. In the case of J1614-2230, such an evolutionary argument of iTauris et al.l (|201l[ ) places 



Table 2. Evolutionary stages for PSR J1614-2230 and PSR J0751±1807 systems. 



ZAMS 


1st RG 


CE SNII 


2nd RG 


now 


J1614-2230 primary 25 M 
secondary 4.5 Mq 
Povb 1 y 


5 x 10 b y 


7 Mq 1.9 Mq 
4.5 Mq 4.5 Mq 

4d 2d 


5 x 10 7 y 


1.97 Mq 
0.50 Mq 

8.7 d 


J0751±1807 primary 15 Mq 
secondary 1.6 Mq 
Porb 1 y 


10 Y y 


5Mq 1.1 Mq 
1.6 Mq 1.6 Mq 

Id 1.5 d 


10 9 y 


1.26 Mq 
0.12 Mq 

6 h 



Formation of millisecond pulsars via accretion with magnetic field 
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a constraint on the accretion time r: r ~ 100 Myr, which results in a pre-accretion mass of 
this NS very close to the observed mass: > 1.9 Mq, and a way to estimate M needed for the 
spin- up (see second term in Eg. 12.11 small r requires large M close to the Eddington limit). 

From Fig. [1] and [5] one can deduce the minimal amount of accreted mass needed to reach 
the observed ~ 300 Hz spin frequency: AA/ min ~ 0.04 ~ 0.05 Mq for r < 1 Gyr. Hence, M in 
such a case is 

M = ^± ~ 10- 9 M /yr t 50 ~\ (3.1) 

T 

where T50 = r/50 Myr. For PSR J0751+1807, our simulations indicate an initial mass of 
~ 1.2M Q ; its value is a function the accretion time r, M and initial magnetic field (see right 
panel in Fig. [2]). 

The influence of the dense-matter EOS can be summarized in terms of its stiffness: the 
stiffer the EOS, the larger the moment of inertia I for a given mass (/ ~ MR 2 ); it requires 
therefore more accretion of angular momentum to spin-up a star to a given frequency in a 
given amount of time (see Eg. 12. ip . While the DH and APR EOSs give very similar results for 
observed frequencies, the BM EOS differs from them since it must be much stiffer (/ is 20% 
higher than for DH EOS) to sustain the softening by hyperons (see e.g., left panel on Fig. [2]). 




M. 




M [M Q 



Figure 1. Currently most massive millisecond pulsar PSR J1614-2230. Left panel: Accreting time 
needed to spin-up the star to presently observed values as a function of accreted mass for DH and 
BM EOSs (Bi = 10 12 G, dotted and dashed lines correspond to lower and upper mass limit, 1.93 
and 2.01 Mq). Right panel: selected spin-up tracks parametrized by different initial masses, magnetic 
fields, accretion rates and times. Lines cross at the current values of the pulsar parameters (BM EOS). 



4. Summary 

We present a versatile method of studying the formation and evolution of millisecond pul- 
sars in the recycling phase. This prescription allows for testing the magnetic field decay and 
magnetic torque models, in the framework of a relativistic thin accretion disk spin-up scenario. 
Coupled with evolutionary data, the information relevant to the SN collapse simulations, spin 
and initial mass distribution of pulsar populations can be drawn from the timing data and 
mass measurements. We point out that from the dense-matter EOS point of view, both the 
low and the high end of the initial mass spectrum is very important, as both the minimal 
and the maximal initial NS mass contains the information about the composition of the mat- 
ter and processes occurring at its birth. Moreover, the disk-magnetic field interaction in our 
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Figure 2. As in Figure 1, but for a low-mass millisecond pulsar PSR J0751+1807. 

model changes dramatically the character of the evolution of the inte rior and the spin of NS, 
as compared to the results where B is neglected ( Bejger et al.l 20l"TbT) . 
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